Entanglement Dynamics in Three Qubit X-States 
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I explore the entanglement dynamics of a three qubit system in an initial X-state undergoing 
decoherence including the possible exhibition of entanglement sudden death (ESD). To quantify 
entanglement I utilize negativity measures and make use of appropriate entanglement witnesses. 
The negativity results are then extended to X-states with an arbitraty number of qubits. I also 
demonstrate non-standard behavior of the tri-partite negativity entanglement metric, its sudden 
appearance after some amount of decoherence followed quickly by its disappearance. Finally, I solve 
for a lower bound on the three qubit X-state concurrence, demonstrate when this bound goes to 
zero, and outline simplifcations for the calculation of higher order X-state concurrences. 
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I. INTRODUCTION 

Entanglement is a quantum mechanical phenomenon 
in which quantum systems exhibit correlations above and 
beyond what is classically possible. As such it is a crucial 
resource for many aspects of quantum information pro- 
cessing including quantum computation, quantum cryp- 
tography and communications, and quantum metrology 
Due to its fundamental, and increasingly practical, 
importance there is a growing body of literature dedi- 
cated to studies of entanglement. Nevertheless, many as- 
pects of entanglement, especially multi-partite entangle- 
ment and its evolution, are in need of further exploration 

i- 

The unavoidable degradation of entanglement due 
to decoherence has severely hampered experimental at- 
tempts to realize quantum information protocols. Deco- 
herence is a result of unwanted interactions between the 
system of interest and its environment. Highly entangled, 
and thus highly non-classical, states may be severely cor- 
rupted by decoherence Q . This is especially troubling as 
these states tend to be the most potentially useful for 
quantum information protocols. An extreme manifesta- 
tion of the detrimental effects of decoherence on entan- 
glement is entanglement suddent death (ESD): in which 
decoherence causes a complete loss of entanglement in a 
finite time d, [f| despite the fact that the system coher- 
ence goes to zero only asymptotically. Much has been 
written about this aspect of entanglement for bi-partite 
systems and there have been several initial experimental 
studies of this phenomenon Q. Fewer studies look at 
ESD, and specifically ESD with respect to multi-partite 
entanglement, in multi-partite systems 04H1- 

A class of two qubit states that are generally known 
to exhibit ESD are the so called X-states [l2[, so named 
due to the pattern of non-zero density matrix elements. 
These states play an important role in a number of phys- 
ical systems and allow for easy calculation of cer- 
tain entanglement measures. In this paper, I explore the 
entanglement dynamics of three qubit X-states in de- 
phasing and depolarizing environments as a function of 
decoherence strength. Previous studies of three qubit X- 



shaped states utilize more restrictive sets of states: GHZ- 
diagonal states [13] and generalized GHZ-diagonal states 
[l5j . Other papers have examined specific examples of 
three qubit X-state entanglement including the effects 
of dephasing on a three-qubit quantum error correction 
protocol [lq ]. 

To quantify entanglement within the three qubit sys- 
tems I utilize the negativity, Nj, defined as the most 
negative eigenvalue of the partial transpose of the density 
matrix jl7| with respect to qubit j. This provides three 
distinct entanglement measures. As a pure tri-partite 
entanglement metric for mixed states I will use the tri- 
partite negativity, which is simply the third root of 
the product of the negativities with respect to each of 
three qubits [3, N (3) = (X1X2X3) 1 / 3 . A mixed state 
with non-zero is distillablc to a GHZ state. It is 

important to note the existence of bound entanglement 
which may be present even if all negativity measures in a 
three qubit system are equal to zero. Thus, when I refer 
to ESD of a state with respect to given negativity metrics 
this should not be confused with separability of the state. 
Nevertheless, besides general interest in the behavior of 
these entanglement metrics, the disappearance of nega- 
tivity plays an important role in quantum information 
protocols in that it indicates that the entanglement of 
the state is not distillable. 

The physical significance of X-states mentioned above 
demands and efficient means of experimentally determin- 
ing the presence of entanglement. This can be accom- 
plished via 'entanglement witnesses.' Three qubit states 
can be separated into four broad categories: separable (in 
all three qubits), biscparable, and there exist two types 
of locally inequivalent tri-partite entanglement (GHZ and 
W-type) [lj|. Reference [2(| provides a similar classifi- 
cation schemes for mixed states each of which includes 
within it the previous classes. These are separable (S) 
states, bi-separable (B) states, W states, and GHZ states, 
which encompasses the complete set of three qubit states. 

Entanglement witnesses are used to determine in which 
class a given state belongs. These observables give a posi- 
tive or zero expectation value for all states of a given class 
and negative expectation values for at least one state 
in a higher (i.e. more inclusive) class. I will make use 
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of specific entanglement witnesses [20J that will identify 
whether a state is in the GHZ\W class (i.e. a state in 
the GHZ class but not in the W class), in which case the 
state has experimentally observable GHZ-type tri-partite 
entanglement. Though the use of entanglement witness 
cannot guarantee that entanglement is not present, it 
does give experimental bounds on whether the entangle- 
ment can be observed. 

The results presented in this paper are (i) the analyt- 
ical determination of various negativity measures for X- 
states of an arbitrary number of qubits including how the 
negativity evolves under decoherence, (ii) the demonstra- 
tion that negativity disappears in finite time for A-states 
subject to different types of decoherence and the (ana- 
lytical and numerical) determination of the decoherence 
strength when this occurs, (hi) the analytical calculation 
of the expectation value of X-states undergoing decoher- 
ence with respect to appropriate entanglement witnesses, 
(iv) the demonstration of the sudden appearance, only at 
non-zero decoherence strength, in some A-states of the 
tri-partite negativity, and (v) the calculation of a bound 
on the three qubit concurrence for A-states and the de- 
scription of how this can be extended to more qubits. In 
addition, I prove in the Appendix that the set of gen- 
eralized GHZ-diagonal states do not cover all possible 
A-states. 



II. THREE QUBIT X-STATES 

There are a number of classes of three qubit states 
whose entanglement properties have been studied and 
whose non-zero density matrix elements form an X 
shape: 



The density matrix of this state is /Ojf(|a| 
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where j = I, ...,4. The most basic is a pure three qubit 
GHZ state with wavefunction \xjj^) = -^(|fc)±|fc)), where 
k is a three bit binary number between zero and seven 
and k is the result of flipping each bit of k. The den- 
sity matrix of this state is px(l/2j, 1/2^-, ±l/2j). Mixed 
states that are diagonal in the basis of these ei ght states 
form the set of GHZ-diagonal states studied in [14[ . The 



basis states have coefficients y A* f° r all < fc < 3, the 
squares of which sum to one. The density matrix ele- 
ments of these states are thus aj = bj = + A^ and 

c j = A fe - A fc- 

A generalized GHZ state is a non-maximally entangled 

state of the form: 



for j = 1,...,4. An incoherent mixture of generalized 
GHZ states where k now ranges from to 7 (as opposed 
to to 3 used in in 

H3) form a generalized GHZ-diagonal 
state. These states are studied in [l5[ and have the form: 

7 

P = XX K P)) («, P) l+Afc («, PWk («. P) I • 



fe=0 



(3) 

For these states the density matrix elements are as fol- 
lows: 



l/?| 2 (At + A,7) + M 2 (A+ + A-) 
aP*{\i-\- k )+a*P{\+-\-) 



(4) 
(5) 
(6) 



(a,/*)) 



a\k)±f3\k). 



(2) 



for 1 < j < N/2 and k = j — 1. However, as shown in 
the Appendix, generalized GHZ-diagonal states do not 
include all possible A-states. This is due to the restric- 
tion of constant a and /3 for all contributing generalized 
GHZ states. 

In this paper I consider A-states that are completely 
general, limited only by the restriction that the state is 
a proper density matrix. For convenience, I will refer to 
the four density matrix elements aj , bj , Cj and c* of the 
X-state as a GHZ-type state. At most, four GHZ-type 
states contribute to each three qubit A-state. 



III. X-STATE ENTANGLEMENT 

An A-state is a mixed state that can be written as a 
sum of GHZ-type states. When a partial trace is taken 
over any one of the three qubits of an A-state the re- 
sulting two qubit matrix is diagonal. This demonstrates 
that the entanglement of an A-state is cither tri-partite 
or biseparable but not completely separable. Before cal- 
culating any specific entanglement metric and studying 
its decay in a given decohering environment, we note that 
an upper bound on the entanglement decay was derived 
in [15| for a number of different decohering environments. 
Though these bounds were calculated for the more lim- 
ited generalized GHZ-diagonal states they appear to be 
appropriate to the states studied in this work. However, 
these upper bounds go to zero only in the limit of com- 
plete decoherence. Thus, the states never exhibit en- 
tanglement sudden death for any entanglement metric 
and the bounds cannot be used to study the ESD phe- 
nomenon. Below, I explore specific entanglement met- 
rics for which I provide analytical solutions to exactly 
calculate the decoherence strength at which the A-states 
exhibit ESD for the given entanglement metrics. While 
ESD of these metrics cannot guarantee separability of the 
A-state it does provide important information concern- 
ing distillability and the ability to determine the presence 
of entanglement . 

To calculate the negativity of a three qubit A-state 
we take the eigenvalues of the partial transpose of the 
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density matrix with respect to one of the qubits. These 
24 eigenvalues (8 for each possible partial transpose) are 
all of the form: 




for all i,j = 1, ...,4 and i ^ j. From these eigenvalues 
one can see how the negativity detects the entanglement 
of an A-state. Let us first assume a GHZ-type state with 
the only non-zero elements a,j , bj , Cj and c* . The eigenval- 
ues which utilize elements <Xj and bj cannot be negative 
(since aj + b 3 ■ = 1 and C, — for all i ^= j). An addi- 
tional three eigenvalues will be equal to — \cj\, demon- 
strating the entanglement in the system. A-states that 
are sums of two GHZ type states have non-zero elements 
aj , aj , bi , bj , c, , Cj ,c*,c*. Such states will again have neg- 
ative eigenvalues — \a\, —\cj\ and two additional possibly 
negative eigenvalues \{a k + b k - \/{a k - b k ) 2 + 4|q| 2 ) 
where k,i = i,j and k ^ £. As more density matrix el- 
ements of the A-state are filled up the eigenvalues tend 
to have the form of these latter two eigenvalues. 
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FIG. 1: (Color online) Eigenvalues of partially transposed X- 
state density matrix with non-zero elements a\ = jt,a2 = 

j2,bi = |,&2 = Y2' Cl = f> aQ d c 2 = |- When tak- 
ing the partial transpose with respect to the third qubit 
there are two 'negative' eigen values which can be c alcu- 
lated from Eq. d | - -L ^625 - 576p(3 + p(p - 3) and 
5j (8 - A/l3-9p(3 + p(p-3))) . Using Eq. [10] we can see that 
only the first of these eigenvalues can be negative and this for 
p < .366. Thus, the state exhibits ESD with respect to N3. 
These two eigenvalues are shown in the plot (solid and dashed 
line respectively). In addition, the two negative eigenvalues 
that constitute iVi and N3 are shown (dotted lines) neither of 
these go to zero in finite time. Thus, the state exhibits ESD 
with respect to the tri-partite negativity at the same value as 
exhibited for N3. 



A. Dephasing Environment 

We now look at the entanglement evolution of the three 
qubit A-states with no interaction between the qubits, 
in an independent qubit dephasing environment noting 
the exhibition of ESD with respect to the negativity. 
The independent qubit dephasing environment is fully 
described by the Kraus operators 

where the dephasing parameter p can also be written 
in a time-dependent fashion, p = 1 — exp(— nt). When 
all three qubits undergo dephasing we have eight Kraus 
operators each of the form Ai = (Ki ® Kj (g> K k ) where 
I = 1,2, ...,8 and i,j,k= 1,2. 

The effect of a dephasing environment of strength p on 
an X-state is to reduce the anti-diagonal elements of the 
density matrix by a factor (1 — p) 3 ^ 2 while leaving the 
diagonal elements constant. To calculate the negativity 
we look at the eigenvalues of the X-state after taking the 
partial transpose with respect to the desired subsystem. 
The relevant eigenvalues are now of the form: 

\ U + bj y / (a,-6 J )2+4|c i P(l-^ . ( 9 ) 
The eigenvalues go to zero when: 



up of single GHZ-typc states (with non-zero elements 
aj + bj i = 1, and Cj): the negativity with respect to any 
one qubit, and thus the tri-partite negativity as well, is 
simply — |cj|(l — p) 3 / 2 . When the X-state is a mixture of 
GHZ type states ESD may be exhibited. Fig. [1] shows a 
sample A-state that is the sum of two GHZ-type states 
that exhibits ESD with respect to the negativity of the 
third qubit, A3. However, the state docs not exhibit 
ESD with respect to N± and N2, they are negative for 
any value of p. The state thus exhibits ESD with re- 
spect to the tri-partite negativity at the same dephasing 
strength as A3. For stronger dephasing no tri-partite 
entanglement is detected. 

The above can be compared to the experimental detec- 
tion capabilities of entanglement witnesses. Appropriate 
entanglement witnesses for A-states are of the sort: 

W k = ^1 - \GHZ{k))(GHZ{k)\ (11) 

where \GHZ(k)) = ^(\k) + |fc)). For an A-state con- 
sisting of a single GHZ-type state the entanglement wit- 
ness W k gives Tr[W kP ] = ±(1 - 4(1 - p) 3 / 2 |c fc |). Thus, 
Wfc loses its ability to detect entanglement at dephasing 
strength p = 1 — 2i / 3 ^ 2/3 ■ The maximum occurs for 
l c fc| = 1/2 in which case the entanglement is no longer 
detected at p = 1— For general A-states the entan- 
glement witnesses give the following: 



(10) 



Based on the above, it is easy to see that ESD with re- 
spect to negativity is not exhibited by A-states made 



Tr[W £ p] = -(3(oj + i>i + aj + bj+a k +b k ) 

+ a e + b e -4(l~p) 3 / 2 \c e \). (12) 



This can be solved for the exact value of p at which the 



4 



entanglement is no longer detected: 
(3(aj + bi + c 



p=l 



b 3 + a k + b k ) + a e + b t ) 2 / 3 



2 4 / 3 |q| 2 / 3 



(13) 

I note that which of the above witnesses is most sen- 
sitive may depend on the intial state and the decoher- 
ence strength and can be determined via a minimization 
process. What is important is that the witnesses de- 
tect purely tri-partite entanglement that does not include 
biseparable but not completely separable entanglement. 



B. Depolarizing Environment 

I now look at an independent qubit depolarizing envi- 
ronment and, as above, explore the entanglement evolu- 
tion of the three qubit X-states. The Kraus operators 
for this environment are: 



Ai 



1 £-1) K w = ~^~ a u 



(14) 



where a w are the Pauli spin operators, w = x,y, z and p 
is now the depolarizing strength. The depolarizing envi- 
ronment affects both the anti-diagonal and diagonal ele- 
ments of the density matrix. The anti-diagonal elements 
are simply reduced by a factor of (1 — p) 3 . The diagonal 
element a* becomes: 

3p 3p 2 p 3 

a, = a.At 

2 4 8 ; 

2 3 



p 



bi — + (at + bj 



(15) 



where if (i,£) — (1,4), (j,k) = (2,3) and vice versa. For 
the b[ elements simply replace each term bi with ai and 
each ai with 67, for I = 1, ...,4. Since the decohering en- 
vironment preserves the X shape of the density matrix 
the eigenvalues of the partially transposed density ma- 
trix follow Eq. [7] and the critical value of p for which 
a given eigenvalue goes from negative to positive can be 
analytically determined. 

When the initial density matrix is composed of only 
one GHZ-type state eigenvalues of the partially trans- 
posed density matrix are the same for each qubit and 
the state exhibits ESD at the same depolarizing strength 
for all Nj and A^ 3 '. When the X-state density matrix 
is a mixture of multiple GHZ-type states ESD may be 
exhibited with respect to specific negativity measures at 
different depolarizing strengths. Figure [2] shows the low- 
est eigenvalue of the partially transposed density matrix 
with respect to each of the three qubits for a sample 
X-state composed of a mixture of GHZ-type states as a 
function of decoherence strength. One eigenvalue is al- 
ways positive (i.e. indicating zero negativity) and two of 
the eigenvalues cross zero (i.e. the state undergoes ESD 
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FIG. 2: (Color online) Eigenvalues of partially transposed 
X-state density matrix with non-zero elements a\ = \ , ai 

1 »- " 1 b s - 3 



g,Gt3 — -§?&4 
To i ^2 q i C'3 



= i, and c 4 = ±. 
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3 = fg , 04 = f , Cl = 

12 , ^ — g , — xq , ^« — 26 ■ The most negative eigen- 
value is shown for the partial transpose of the initial X-state 
density matrix taken with respect to the first (dotted line), 
second (solid line), and third (dashed line) qubit. The state 
exhibits ESD with respect to the latter two negativities at 
dephasing values p ~ .0585 and .0317, respectively. The in- 
set shows the tri-partite negativity. Note that this tri-partite 
entanglement measure goes from positive (indicating the lack 
of entanglement) to negative (indicating the presence of 
GHZ distillable entanglement), back to positive. This entan- 
glement behavior arises from the fact that the single qubit 
negativities cross zero at different decoherence strengths. 



with respect to the single qubit negativities) at different 
decoherence strengths. For low values of p two of the low- 
est eigenvalues are negative demonstrating the presence 
of entanglement. However, there is no measurable GHZ 
distillable tri-partite entanglement as measured by N^ 3 \ 
For slightly higher values of p there is a small region for 
which only one of the eigenvalues is negative. Now 
becomes negative showing a sudden birth of (GHZ distil- 
lable) tri-partite negativity. As p increases further, the 
state exhibits ESD with respect to all single qubit nega- 
tivities and 7V( 3 ) (and A^) becomes positive. This sort of 
A^ 3 ' behavior indicates that there is only a small region 
of decoherence strengths (which does not include p = 0) 
for which we can be sure there exists GHZ-distillable en- 
tanglement. Such behavior, going from positive to neg- 
ative and back, cannot occur when the X-state is com- 
posed of only two GHZ-type states. This is because two 
of the single qubit negativites are equal, the partial trace 
with respect to two of the qubits give the same set of 
eigenvalues. The sign of iV^ 3 ) is thus determined solely 
by the eigenvalues of the partially transposed state with 
respect to the third qubit. An example of a state exhibit- 
ing the sudden birth of A^ 3 ) entanglement followed by an 
exhibition of ESD is portrayed in the inset of Fig. [5] 

To test for the presence of purely tri-partite entangle- 
ment in the depolarizing system wc look at the expec- 
tation value of the state with an entanglement witness. 
Using the witnesses defined above we note that witness 
Wj for a state depolarized with a strength p gives: 



Tr[W jP ] = \(p 2 -2p + eX^Oi + hi) 



-(3p 2 - 6p - 2) (cy +bj) + (p- l) 3 c/.16) 
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This equation can then be solved for the critical decoh- 
ernce strength at which entanglement will no longer be 
detected. 



IV. X-STATES WITH MORE QUBITS 

The particular matrix structure of the A-state allows 
us to extend our results beyond three qubits. A matrix 
with non-zero elements in an X shape can be block diago- 
nalized with blocks of size 2x2. Assuming the A-matrix 
elements along the diagonal are di,...,djv, and the ele- 
ments along the anti-diagonal are ei, ...,ejy (starting at 
the top right), the mth 2x2 block along the diagonal 
has elements: 



dm &m 
CN-m+l dN-m+1 



(17) 



where 1 < m < N/2. Thus, the eigenvalues of any di- 
mension A-shaped matrices are simply the eigenvalues of 
the 2x2 blocks which are 



the lower bound goes to zero. Once the lower bound docs 
go to zero, there is no longer a guarantee that entangle- 
ment is present. 

For initial density matrices that are A-states only six 
of the 18 contributing terms to the three qubit con- 
currence are non-zero. More specifically, only the two 
SO (4) generators with elements on the anti-diagonal 
contribute to each of the three bi-partitc concurrences. 
I will refer to the SO (4) generator with anti-diagonal 
( — 1,0,0,1) as , and the SO(4) generator with anti- 
diagonal (0, —1, 1, 0) as L l 2 . For AT-states the four eigen- 
values that make up each of the six terms are of the form: 



A 



ij\k 



±2^ 



(21) 



for two different values of m. For the partition 12 1 3 
and SO (4) generator 1 = 1 the contributing terms have 
r?i = 1,2. For the generator £ = 2,m = 3,4. Simi- 
larly, for the partition 23 1 1 we find I = l,m = 2,3 
and I = 2,m = 1,4. Finally, for the 13 1 2 partition 
I = l,?n = 2,4 and £ = 2,777 = 1,3. Given these eigen- 
values the three-concurrence can be easily computed. 



-^{d m + dN-m+l ± y/ {dm — ^AT-m+l) 2 + 4e r „ejV-m+l ) ■ 

(18) 

When calculating the negativity a partial transpose of 
the density matrix must be taken. This has the effect 
of rearranging only the elements along the anti-diagonal 
while preserving the A shape. Thus, the eigenvalues of 
the partial transpose of an A-state of any dimension have 
the form of Eq. [7] and the negativity is easily calculated. 



V. THREE-QUBIT CONCURRENCE 

In this section I derive an explicit expression for a lower 
bound of three qubit mixed state concurrence as defined 
in [2l| for A-states. A general expression for a lower 
bound on the three-qubit concurrence is: 



73 



\ 1=1 



Each of the three bi-partite concurrence terms C^' k is 
given as the sum of the six terms: 



c{0, JX 



(20) 



where \\ arc the non-zero eigenvalues of p = 

pSg p*S % l in descending order. The operators 

are given by Sf^ k = Lif <g) L§ where L 1 ^ is one of six 
generators of the group SO (4) operating on qubits i,j, 
and Lq is the generator of SO (2), the Pauli matrix a y , 
operating on qubit k. This lower bound on mixed state 
concurrence has been calculated for some simple A-states 
in [22! . Here I look to extend these results and note where 



A. Dephasing Environment 

For A-states composed of only one GHZ-type state an 
exact calculation for T3 in a dephasing environment yields 
the maximum between and: 



a t -af + Ci(-uii + 27/ ) - \J a, - af + c^-Wj - 27? ) 

(22) 



where, 



= Ci(p - l) 3 
ji = ai(a,i - l)(p - if. 



(23) 



Eq. [22] goes to zero only in the limit of p — » 1. Thus, 
the lower bound cannot go to zero for a GHZ-type state, 
some entanglement will always be present. In fact, the 
lower bound cannot go to zero unless the A-state is com- 
posed of four GHZ-type states. This is because T3 is a 
summation of terms and can go to zero only if each term 
goes to zero. Each one of these (six) terms consists of four 
eigenvalues, two for each of two m values. If the eigen- 
values of one of the 777 values are zero (which will happen 
if Ci and or bi = 0) the term will have the form of Eq. 
1221 Thus, that term, if not initially zero, will remain non- 
zero until p = 1. This behavior is demonstrated in Fig. 
[3] and should be contrasted with the negativity and tri- 
partite negativity measures. The negativity of X-states 
can go to zero in a dephasing channel when the A'-state 
is composed of only two GHZ-type states. The reason for 
this is that the negativity can be defined with respect to 
only one of the qubits (for example A3) which may go to 
zero while negativity measures with respect to the other 
qubits do not. The three qubit concurrence, however, is a 
sum over all terms and therefore cannot go to zero unless 
each bi-partitc concurrence term goes to zero. 
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FIG. 3: (Color online) Top row: Three qubit concurrence for 
X-states composed of one GHZ-type state with a» = bi = 
1/2. Left: in a dephasing environment T3 never goes to zero. 
Right: in a depolarizing environment T3 goes to zero when 
\{p(p — 2) + 8ci(l — p) 3 ) = 0. Bottom: example of three-qubit 
concurrence evolution as a function of dephasing strength. 
The chosen X-state density matrix is: ai = | — e, 0,2 = ^,03 = 
e,a 4 = Teibi = ^,62 = ^,63 = j^,b4 = g,ci = ^,c 2 = 
§,C3 = |, and C4 = T3 is shown for e = (solid line), 
e = 1/128 (dotted line), e = 1/32 (dashed line), and e = 1/16 
(chained line). Note that when G13 = T3 does not go to zero. 
Adding any finite amount to 03 allows T3 to go to zero. 

B. Depolarizing Environment 

The bound on the three-qubit concurrence for X-states 
composed of one GHZ-typc state in a depolarizing envi- 
ronment gives the maximum between and: 

*_.L_2i_ v^h (24) 

4 V 64 4 V 64 4 V ' 
where, 

q t = (p - 2) P y / 4(p - l) 2 (a> - a?) - p(p - 2) (25) 

7.I = (a i (p-2) 3 + (a. 1 -l)F 3 )((a l -l)(p-2) 3 + aip 3 ). 

T3 for this state is shown in Fig. [3] for an initial state 
ai = bi = 1/2. For the depolarizing environment T3 can 
go to zero even for X-states composed of only one GHZ- 
type state. 

The matrix p for a X-state density matrix with any of 
the six SO (4) generators retains its X shape, having four 
diagonal and four anti-diagonal elements (and thus have 
only four non-zero eigenvalues as noted in [2l|). The 

X is 'balanced' when Sj 7 ''* is anti-diagonal. I use the 
term 'balanced' as follows: any diagonal X-matrix ele- 
ment d m that is non-zero has a non-zero counterpart el- 
ement e TO , where I have used the notation of Eq. 1171 An 
unbalanced X-matrix will have non-zero elements whose 
counterparts are zero. When a balanced X-matrix is 
block diagonalized non-zero 2x2 blocks will have four 



non-zero elements leading to non-degenerate eigenvalues 
like those of Eq. [T51 Block diagonalized unbalanced X- 
matriccs will have a zero in one of the off-diagonal ele- 
ments of the 2x2 diagonal block as can be noted from 
Eq. [TTl In the unbalanced case the eigenvalues are then 
simply the diagonal elements of the block which are of 
the form aibj,aibj,akbi,akbi. These elements (eigenval- 
ues) are degenerate and thus these bi-partite concurrence 
terms equal zero. 

C. n-Qubit Concurrence 

As mentioned above, only the two SO (4) generators 
with elements on the anti-diagonal, contribute to the 
three qubit concurrence. This is because only these two 
generators have anti-diagonal elements which lead to bal- 
anced X-matrices p. The other generators lead to unbal- 
anced p matrices whose eigenvalues are simply its diago- 
nal elements. The eigenvalues are each doubly degenerate 
meaning that these terms will not contribute to the con- 
currence. The above allows us to simplify calculations for 
higher qubit concurrences of X-states. The only terms 
necessary to calculate are those that utilize anti-diagonal 
S/" " matrices. Thus, for four qubits there would be 
four terms from the SO (4) <g> SO (4) generators for each 
of the three balanced partitions (two qubits on each side 
of the partition) and an additional four terms from the 
SO(8)®SO(2) generators for each of the four unbalanced 
partitions (partitions of three and one qubit). This gives 
a total of 28 terms which should significantly simplify 
these calculations. 



VI. CONCLUSIONS 

In this paper I have studied the entanglement dynam- 
ics for three qubit X-states in both dephasing and de- 
polarizing environments. To do this I have analytically 
calculated the eigenvalues of partial transposes of the X- 
states which allows for easy determination of the negativ- 
ity. Since the dephasing and depolarizing environments 
retain the density matrix X shape one can calculate 
which initial states will exhibit ESD with respect to the 
negativity measures and at what decoherence strength. I 
noted that the tri-partite negativity, a tri-partite entan- 
glement measure which is sufficient to ensure GHZ dis- 
tillability, can exhibit non-standard behavior for certain 
X-states: its appearance only at non-zero decoherence 
strength followed by its sudden disappearance. In addi- 
tion, I explored the detection capability of entanglement 
witnesses sensitive to tri-partite entanglement. As with 
the negativity, the expectation value of the X-state with 
respect to the entanglement witness can be solved ana- 
lytically and are vital in assessing potential experimental 
studies. These results are extended to systems made of 
arbitrary numbers of qubits. Finally, I analytically solved 
for the relevant terms of a lower bound on the three-qubit 
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concurrence for an A-state, demonstrated when it goes to 
zero in dephasing and depolarizing environments. This 
method may be useful for calculating concurrences for 
larger numbers of qubits. 

It is a pleasure to thank G. Gilbert and S. Pap- 
pas for helpful feedback and acknowledge support from 
the MITRE Innovation Program under MIP grant 
#20MSR053. 



or (3 = 

If either of these is true then all other c m for m ^ j 
must also equal zero. 



A fc = X k and A r = h 



If this is true a 
2(|/3| 2 A fc + H 2 A r ) where A 



i - 2(|a| 2 A fe 



| 2 A^) and bj = 
- XZ . Therefore, 



if in addition bj = (which would require A^ 
Xr = 0), dj must equal zero. 



Appendix A: General X-States 

As mentioned in the main part of the paper, a previ- 
ously studied set of states with an A shaped density ma- 
trix is the generalized GHZ-diagonal states [l5|. In this 
Appendix I prove that generalized GHZ-diagonal states 
do not include all possible X-states by constructing an 
explicit state with an X-shaped density matrix that is 
not part of the aforementioned set. 

Generalized GHZ-diagonal states with n qubits have 
the form 

N-l 

P= E A+|^ + («,/?))(^ + («,/3)|+A-|^T(«,/3))(^7(a,/3)|, 

fe=0 

(Al) 

where N = 2" is the Hilbert space dimension. The den- 
sity matrix elements of these states using the notation of 
Eq. [I] are as follows: 

aj = H 2 (A+ + A-) + |/?| 2 (A+ + A_) ( A2 ) 
b 3 = |/?| 2 (A+ + A,;) + H 2 (A+ + A-) (A3) 
Cj = ap*{Xi-Xl) + a*p{X+-XZ) (A4) 

for 1 < j < N/2 and k= j — 1. 

We now construct an X-state that is not a generalized 
GHZ-diagonal state. Let us set Cj — 0. There are then 
three possible solutions for Eq. A4: 



Let a* (3 ■■ 
tion SC£ = 



4) 



re where r, 9 are real. Then the frac- 
As mentioned in the main part 
of the paper, the coefficients A are all real forcing 
e 2lB to be real and 9 = 0, mn for all integers m. 
Therefore, a* f3 and aj3* must both be purely real 
or purely imaginary. 
We can now explicitly construct a two qubit X-state 
that is not part of the set of generalized GHZ-diagonal 
states by setting c\ = and violating each of the three 
conditions listed above. Such a state can have the form: 



PC 



( 0,1 




V 





a 2 





J4> 



re- l(t > b 2 
/ 



(A5) 



where r, 9 are real and r 2 < a 2 b 2 guarantees the density 
matrix has positive eigenvalues. In addition, pc must 
be trace 1, 01+02 + 62 = 1, and its purity must be 
a 2 + a 2 + 6 2 + 2r 2 < 1. 

In Eq. IA5l ci = 0, yet c 2 7^ 0, indicating that a, (3 ^ 0. 
Furthermore, b\ = while a\ does not, indicating that 
Xq = Xq and A;j~ = A^ cannot both be true. Finally, c 2 = 
re 1 ^ need not be real nor purely imaginary indicating that 
(Aq — Aq~) 7^ (A^ — A3"). Thus, the state pc is not 
part of the set of generalized GHZ-diagonal states though 
it certainly is an A-statc. 
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